A theoretical analysis on highly spin-polarized transport of iron nitride Fe 4 N 
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In order to propose a ferromagnet exhibiting highly spin-polarized transport, we theoretically 
analyzed the spin polarization ratio of the conductivity of the bulk Fe4N with a perovskite type 
structure, in which N is located at the body center position of fcc-Fe. The spin polarization ratio is 
defined by P — (<tj — <7j)/(oj +0|), with cr-m) being the conductivity at zero temperature of the up 
spin (down spin). The conductivity is obtained by using the Kubo formula and the Slater- Foster 
tight binding model, where parameters are determined from the least-square fitting of the dispersion 
curves by the tight binding model to those by the first principles calculation. In the vicinity of the 
Fermi energy, \P\ takes almost 1.0, indicating perfectly spin-polarized transport. In addition, by 
comparing Fe4N to fcc-Fe (Fe4No) in the ferromagnetic state with the equilibrium lattice constant 
of Fe4N, it is shown that the non-magnetic atom N plays an important role in increasing \P\. 

PACS numbers: 72.25.Ba 
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Recently, highly efficient spin-electronics devices op- 
erating at room temperature have been extensively de- 
veloped for applications to the magnetic memory and 
the magnetic sensor. A typical device has ferromagnetic 
tunnel junctions consisting of a ferromagnetic electrode 
(FME)/insulator/FME, which exhibits a large magne- 
toresistance (MR) effect. 1 ^ 3 ! 4 ! 5 ! 6 ! 7 The efficiency of the 
MR effect is often defined by the MR ratio = (R AP — 
R P )/R Ar , with i? P and R Ar being the resistance of 
the parallel and anti-parallel magnetization configura- 
tions of FMEs, respectively. Experimentally, regarding 
junctions with an electrode of a half-metallic ferromag- 
net, C02Cro.6Feo.4Al/Al-O/CoFe/NiFe/IrMn/Ta junc- 
tions exhibited an MR ratio of 16 % at room temper- 
ature^ and Co 2 MnSi/Al-0/Co 75 Fe25 had an MR ra- 
tio of 70 % at room temperature^ Regarding junc- 
tions of electrodes with usual ferromagnets, single-crystal 
Fe(001)/MgO(00T)/Fe(001) junctions exhibited MR ra- 
tios of 88 %2i and 180 %4 at room temperature. Fur- 
thermore, CoFeB/MgO/CoFeB junctions achieved MR 
ratios of 260 %& and 355 %£ (the world's highest value) 
at room temperature, although the crystal structure of 
CoFeB and the role of the light element B on the spin- 
polarized transport have not been clarified yet. Gen- 
erally, the MR ratio becomes large with increasing the 
spin polarization of the conduction electron in the FME. 
In the future, an FME with more highly spin-polarized 
electrons, which would result in larger MR ratios, will be 
strongly desired from the viewpoint of the development 
of highly efficient MR devices. 

Towards proposal of such an FME, we extracted an 
idea to obtain electrodes exhibiting the highly spin- 
polarized transport at room temperature. We found 
that ferromagnets consisting of magnetic elements and 



light elements, such as CoFeB, might be very useful as 
the electrodes. 

We consider Fe as a representative magnetic element 
in this idea. We define the spin polarization (SP) ratio 
about the density of states (DOS) of the bulk system, 
P oos (E), as 

D,(E)- Dl (E) 

D^(E)+D l (E)' W 

with D S (E) being the DOS of spin s (=f or j) for the 
energy E. Then, |P DOS (£ P )| at the Fermi energy E F of 
fcc-Fe in the ferromagnetic state is evaluated to be about 
0.7 [see Fig. 2(c)], which is about 2.3 times as large as 
that of bcc-Fe^ where the most stable state of bcc-Fe 
(fcc-Fe) is ferromagnetic (not ferromagnetic^) . This may 
indicate that the highly spin-polarized transport is re- 
alized for materials closely related to fcc-Fe. The con- 
ductivity, including the velocity of electrons, will give an 
answer to the realizability of such transport. 

We focus on a ferromagnet containing fcc-Fe and 
a light element, Fe4N with a perovskite-type struc- 
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in which N is located at the body cen- 
ter position of fcc-Fe. This ferromagnet has a Curie tem- 
perature of 761 K>i£ It should be noted that studies on 
the spin-polarized transport of Fe4N have scarcely been 
performed so far, although other properties have been ex- 
perimentalljiiSiiiii^ and theoretically 13 - 14 - 15 investigated. 
We are particularly interested in not only the SP ratio 
on the transport of Fe 4 N but also the role of the light (or 
non-magnetic) element N on the transport. 

In this paper, we analyzed the SP ratio of conduc- 
tivity of Fe4N in order to elucidate the spin-polarized 
transport. The conductivity was obtained for each spin 
and each orbital of the bulk system using the first princi- 
ples (FP) calculation and the tight binding (TB) model 
calculation. Consequently, we found that Fe4N exhibits 
an extremely highly spin-polarized transport and that N 
plays an important role in the transport. 
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A calculation method is introduced to obtain the con- 
ductivity and the SP ratio. The method is a combination 
of (i) the FP calculation^ (ii) the TB modelj±£ and (iii) 
the Kubo formula>i^ Details of each are given below: 

(i) The FP calculations are performed by the Vi- 
enna Ab-initio Simulation Package (VASP) codei& 
based on the spin-polarized density functional the- 
ory, where we employ the generalized gradient ap- 
proximation of Perdew and Wang (PW91) 19 and 
ultrasoft pseudopotentials to describe the core elec- 
trons. The cutoff energy for the plane wave basis is 
237.51 eV for Fe and 348.10 eV for Fe 4 N, and the 
Monkhorst-Pack sel3 of 8x8x8 k points is used. 

(ii) The Slater-Koster TB modelii is used with tak- 
ing into account the 3d, 4s, and 4p orbitals for Fe 
and the 2s and 2p orbitals for N as well as interac- 
tions up to the third-nearest neighbor atoms. The 
Slater-Koster parameters of potential energies and 
transfer integrals^ are here determined from the 
least-square fitting of the dispersion curves by the 
TB model to those by the FP calculation at the 
equilibrium lattice constant. The fitting is done in 
energy regions from the lowest energy to E F + 5 eV 
for bcc-Fe and fcc-Fe and from the lowest energy to 
E F + 2 eV for Fe 4 N. The number of parameters is 
34 for bcc-Fe and fcc-Fe and 87 for Fe 4 N. 

(iii) The Kubo formula^ and the Slater-Koster TB 
model with the determined parameters are used 
to calculate the conductivity at zero temperature. 
In this calculation, we utilize the theory given by 
Tsymbal et alS^ The total conductivity of spin s 
(=T or I), o- s (E), is written by, 
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where <Ji :S (E) is the conductivity of orbital i (=4s, 
4p, 3d orbitals, and so on) and spin s. This <Ji^ s (E) 
includes the velocity of electrons and the Green's 
function. The Green's function has a single param- 
eter 7 in the second-order self energy due to the 
weak electron-impurity interaction, where 7 char- 
acterizes the degree of electron-impurity scatter- 
ing. In detail, 7 is related to the lifetime of the 
electron of orbital i and spin s, -77^, via l/V,- lS = 
(27r/h)j 2 D i}S (E), with Di )S (E) being the partial 
DOS for i and siSi At present, 7 is set to be 0.5 
eV, which was previously chosen so as to reproduce 
the resistivity of copper^ We also use the diago- 
nal approximation for the self energy, in which the 
scattering of electrons due to impurities is allowed 
for the same energy levels but forbidden between 
different energy levelsi^i The SP ratio is then de- 
fined by 



P(E) = 



^(E)-a l (E) 
a r (E)+a L (E)' 
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We first compare the results of bcc-Fe by the present 
method with the previous ones. The equilibrium lattice 
constant obtained by using the FP calculation is esti- 
mated to be 2.84 A with an error of about 1 % for the 
experimental value of 2.87 Ami The total conductivity at 
E F of the down spin ai(E F ) is larger than a^(E F ). The 
SP ratio at E F , P(E F ), therefore takes a negative value, 
—0.20, and it qualitatively agrees with the previous re- 
sult of about —0.26, 21 which was obtained by a similar 
method combined with the FP calculation within the lo- 
cal density approximation. 

In the following, we investigate the equilibrium lattice 
constant, D S (E), and P uos (E) for Fe 4 N using the FP cal- 
culation. The equilibrium lattice constant is evaluated as 
3.810 A, which has an error of less than 1 % for an ex- 
perimental value of 3.795 A^ As shown in Figs. 1(a) 
and 1(c), D L (E F ) is higher than D r (E F ). Partial DOSsM 
schematically illustrated in Fig. 1(b) show that 3d or- 
bitals are dominant around E Fl and 4s and 4p (4s-4p) 
orbitals are mainly located in an energy region higher 
than E Fl while each orbital of N atoms mostly exists in 
an energy region lower than E F . Furthermore, P DOS (E F ) 
obtained by using D S (E F ) of the FP calculation is eval- 
uated to be —0.6 [see Fig. 1(c)]. 

Using the TB model with the parameters determined 
from the fitting of the dispersion curves, we obtain D S (E) 
and P uos (E) for — 1 eV < E < 1 eV. As seen from Fig. 
1(c), D S (E F ) and P uos (E F ) of the TB model agree well 
with the respective ones of the FP calculation. 

With the use of the TB model and the Kubo formula 
we calculate a l , s (E), <t s (E), and P(E) for — 1 eV < 
E — E F < 1 eV. The results are shown in Fig. 1(d). 
For E < E F , a 3d ^(E) becomes relatively large owing 
to the high DOS of the 3d orbitals of the up spin. 
For E - E F > 0.6 eV, the 4s-4p orbitals of the up 
spin contribute strongly to <7f (E) in spite of their low 
DOSs because their orbitals have large velocities. For 
eV < E — E F < 0.6 eV, each ai^(E) has a pronounced 
valley reflecting the low DOS of the up spin. In this en- 
ergy region, although the DOS of the up spin is actually 
lower than that of the FP calculation, the qualitative be- 
havior of ai^(E) appears to be valid because the DOS of 
the FP calculation has very few components of the 4s-4p 
orbitals and it is low. On the other hand, each <7;^(.E) 
is almost flat, and the 3d orbitals of the down spin con- 
tribute largely to &i{E) because of their high DOS. The 
total conductivity at E F of the down spin ai(E F ) is much 
larger than a^(E F ). The SP ratio P(E F ) therefore takes 
almost —1.0, indicating perfectly spin-polarized trans- 
port. The magnitude of the SP ratio |P(-E F )| is about 
5.0 times as large as that of bcc-Fe. 

In order to clarify the effect of an N atom on the trans- 
port, we investigated D S (E), P DOS (E), a ijS (E), a s (E), 
and P{E) of fcc-Fe (Fe 4 No) in the ferromagnetic state 
with an equilibrium lattice constant of Fe 4 N. Regard- 
ing D S (E) around E F obtained using the FP calculation, 
Di(E) is much higher than D^(E) [see Figs. 2(a) and 
2(c)]; for further details, low and broad DOSs of the 4s- 



3 



4p orbitals of the up spin and a high DOS of the 3d 
orbitals of the down spin are observed [see Fig. 2(b)]. 
The SP ratio of DOS, ^dos(-Ef), is evaluated to be about 
—0.7. Figure 2(c) also shows D S (E) and P uos (E) for 
— 1 eV < E < 1 eV; these values are obtained using the 
TB model with the determined parameters. It is found 
that D S {E F ) and P DO s(^f) of the TB model agree fairly 
well with the respective ones of the FP calculation. With 
the use of the TB model, we calculate <Ji yS (E) and a s (E) 
[see Fig. 2(d)]. In spite of the low DOSs of the 4s-4p 
orbitals of the up spin, their orbitals contribute strongly 
to <Tj (E) . This tendency originates from the large veloci- 
ties of the 4s-4p orbitals. Moreover, a 3d ^(E) is relatively 
large reflecting the high DOS of the 3d orbitals of the 
down spin. The total conductivity at E F of the up spin 
a-\(E F ) then becomes larger than ai(E F ) even though 
D^(E F ) is much lower than Di(E F ). The SP ratio P{E F ) 
has a positive sign, opposite to the sign of P DOS (E F ). 
This sign of P(E F ) is also opposite to that of Fe4N. The 
magnitude of the SP ratio \P(E F )\ is evaluated to be 0.4, 
which is between the \P(E F ) \ of Fe4N and that of bcc-Fe. 

On the basis of the above investigations and the evalu- 
ated Slater-Koster parameters, we discuss the role of the 
N atom on the highly spin-polarized transport of Fe4N 
using schematic illustrations of partial DOSs [see Figs. 
1(b) and 2(b)]. Note that by merely introducing the non- 
magnetic atom N to the body center position of fcc-Fe, 
\P(E F )\ becomes about 2.5 times as large as that of fcc- 
Fe and the sign of P(E F ) changes. In the present study, 
we find that, by adding N to fcc-Fe, 4s-4p bands of Fc 
are raised to a higher energy region by a large magnitude 
of transfer integrals between the 4s-4p orbitals of Fe and 
the 2s and 2p orbitals of N, while 3d bands do not change 
significantly owing to the small magnitude of the trans- 
fer integrals between the 3d orbitals of Fe and the 2s and 
2p orbitals of N. These behaviors are explained by con- 
sidering the bonding-antibonding states formed by the 
transfer integrals, which correspond to overlaps between 
both orbitals. A large portion of 4s-4p bands exists in 



an energy region higher than 3d bands, and a small por- 
tion of them is located in the energy region of 3d bands 
because of hybridizations with 3d orbitals, and is then 
spin-polarized there. Hence, comparing Fe^N with fcc- 
Fe, the proportion of 4s-4p bands is small in the vicinity 
of E F , and, in particular, the proportion of the up spin is 
extremely small; the 3d bands of the down spin are dom- 
inant there. Therefore, the magnitude of the SP ratio 
|P(£' F )| increases, and P(E F ) becomes negative. 

Finally, although it would be difficult to derive the 
general properties of phenomena from the present study 
alone, it seems possible that various ferromagnets con- 
sisting of magnetic elements and light elements have a 
large magnitude of SP ratios according to the same mech- 
anism as that of Fe4N. In addition, when such ferromag- 
nets are used as the FMEs, FME/insulator/FME junc- 
tions may exhibit large MR ratios, although the MR ra- 
tios are often influenced by the interfacial states and ma- 
terials of the insulator. In fact, CoFeB electrodes bring 
about a very large MR effect in CoFeB /MgO/ CoFeB 
j unctions 

In conclusion, |P(£' F )| of Fe4N was evaluated to be 
almost 1.0, which was about 5.0 times as large as that 
of bcc-Fe and about 2.5 times as large as that of fcc- 
Fe. In comparison with fcc-Fe, it was shown that the 
large magnitude of the SP ratio originated from the con- 
tribution to the transport of 3d bands, which was en- 
hanced by introducing N. We anticipate that Fe4N will 
become an electrode with a high efficiency of spin in- 
jection. Furthermore, various ferromagnets consisting of 
magnetic elements and light elements may exhibit highly 
spin-polarized transport due to the present mechanism. 
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FIG. 1: Density of states, conductivities, and spin polariza- 
tion ratios of Fe4N: (a) D S (E) calculated using the FP calcu- 
lation, (b) Schematic illustration of partial DOSs. The 4s-4p 
components are partly covered by the 3d ones [see darkish 
parts in 3d components], (c) D 3 (E) and Pdos(E) for —1 eV < 
E - E F < 1 eV of the TB model and the FP calculation, (d) 
(Ti, s (E), a s (E), and P(E) for -1 eV < E - E F < 1 eV. Here, 
C2s,s(E) and a 2p , s (E) for N are not shown because of their 
small values. 
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FIG. 2: The same as Fig. 1 for fcc-Fe (Fe 4 N ) in the ferro- 
magnetic state with the equilibrium lattice constant of Fe4N: 
(a) D s (E) calculated using the FP calculation, (b) Schematic 
illustration of partial DOSs. The 4s-4p components are partly 
covered by the 3d ones [see darkish parts in 3d components], 
(c) D S (E) and P DO s(E) for -1 eV < E - E F < 1 eV of the 
TB model and the FP calculation, (d) a iyS (E), a s {E), and 
P{E) for -1 eV < E - E F < 1 eV. 



